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We present the hrst representation of the general term of the Rayleigh-Schrodinger series for 
quasidegenerate systems. Each term of the series is represented by a tree and there is a straight- 
forward relation between the tree and the analytical expression of the corresponding term. The 
combinatorial and graphical techniques used in the proof of the series expansion allow us to derive 
various resummation formulas of the series. The relation with several combinatorial objects used 
' for special cases (degenerate or non-degenerate systems) is established. 
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I. INTRODUCTION 



Rayleigh-Schrodinger (RS) perturbation theory is a venerable technique to calculate the eigenvalues and eigenvectors 
of H = Hq + V from those of Hq. It was created in 1894 by Lord Rayleigh to describe the vibrations of a string [l[ 
and adapted to quantum mechanics by Schrddinger in 1926 [2]. RS perturbation theory has been used in all fields of 
quantum physics (particle, atomic, molecular, solid-state physics) and quantum chemistry. 

In most textbooks, the RS series deals with the perturbation of a single nondegenerate initial state. However, 
in many practical applications, the initial state is either degenerate or quasidegenerate (i.e. several states are close 
in energy) and the basic RS approach breaks down. Quasidegenerate perturbation theory is widely used to set up 
effective Hamiltonians 3] , for example the spin Hamiltonians of molecular chemistry and solid state physics [3, H| , or 
to deal with the quantum electrodynamics of atoms [f| 0] • 

The terms of the RS series are notoriously complex. Even mathematical physicists of the stature of Reed and Simon 
admit that the terms of the RS series are "quite complicated" and that "the higher order RS coefficients arc hard to 
compute" (see ref. [H, p. 8 and 18). Computer programs are available to build the terms of the RS series [H-flU, but 
their results are intricate expressions that do not exhibit any obvious structure and that can hardly be used to carry 
out resummations of the series. 

In this paper, combinatorial physics is used to provide the first representation of the general term of the Rayleigh- 
Schrodinger (RS) perturbation theory for quasidegenerate systems. Each term is written as a tree that faithfully 
reflects its algebraic structure. In particular, trees illustrate the recursive structure of RS terms, that is used to prove 



H | properties of the RS series. 

The purpose of this paper is to describe the general term of the RS series and to illustrate the power of our 
combinatorial approach by deriving a number of possible resummations of the RS series. Our aim is to set up the 
tools for performing such resummations and not to discuss when some resummations are more convenient than others. 
As a consequence, no numerical example is given. 

When the system is not quasi-degenerate (i.e. when it is either fully degenerate or non-degenerate), several graphical 
representations of the terms of the RS series have been proposed. We give the bijection between these representations 
and our trees. 



II. RAYLEIGH-SCHRODINGER SERIES 



In the most general setting, we consider a model space M spanned by N eigenstates |z) of H , with energies e^. We 
assume that all of the model space are separated from the rest of the spectrum of Hq by a finite gap. The projector 
onto the model space M is P = where i runs over the basis states of M. The wave operator f2 transforms N 

states \4>i) of M into eigenstates \<&i) = ft\4>i) of H: Hfl\(j)i) = Eifl\(j)i). We recently described [3| the way to choose 
the states | </>,). 
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We assume that Pfl = P and flP = fl. Then, the eigenvalues of the effective Hamiltonian H c h = PHfl are the 
eigenvalues £J, of H. In other words, eigenvalues of H can be obtained by diagonalizing the M x M matrix H g. 
Lindgren and Kvasnicka showed [HI EB] that 

[fi,ir ] = vn-nvn. (i) 

The recursive solution of this equation gives a series expansion for SI which is the RS series for the wavefunction. 
However, it is not obvious that eq. ([T]) has a recursive solution. Indeed, we must be able to solve the equation 
[X, Hq] = C for X. In general, the operator Sylvester equation AX — XB = C, where A and B are self-adjoint, 
has a unique solution if and only if A and B have no common eigenvalue pj], Hal- Clearly, this is not the case if 
A = B = Hq, so we recast the equation by defining x = SI — P, where \ = QxP ■> with Q = 1 — P. Thus, eq. (JTJ 
becomes an equation for x- 

[X,H ] = QVP + QVx-xVP-xVx- (2) 

This equation is a matrix Riccati equation, for which various solution methods have been proposed in the literature fl§h 
[2lj |. Solutions exist also for its time-dependent form [H, |23j]. Now, [x, H ] = xHq — Hqx = xPHq — QH x, because 
P and Q commute with Hq . We obtain a Sylvester equation with A = QHq and B = PHq and its solution is unique 
because we assumed that there is a finite gap between the states of the model space and the rest of the spectrum. 



A. Combinatorial analysis 

Because of its importance, the basic RS series has been dealt with in hundreds of papers. As a starting point 
to discover the proper combinatorial structure of the RS series, we enumerate the number of its terms: it has one 
term of order 0, one term of order 1, then 2, 5, 14, 42, 132, 429 terms of order 2, 3, 4, 5, 6 and 7, respectively. 
The number of terms of order n is the Catalan number C n = jjj^ppp- The Catalan sequence is well- know in 
combinatorics: "The Catalan sequence is probably the most frequently encountered sequence that is still obscure 
enough to cause mathematicians [. . . ] to expend inordinate amounts of energy re-discovering formulas that were 



worked out long ago." [24|. An entire book is devoted to Catalan numbers [25]. They enumerate over 200 combinatorial 
objects [2(| and in the physics literature, several of these objects were used to represent the general term of the RS 
series: Dyck paths [27j |. non-crossing partitions (28|, bracketings [HI, [29|, [3(| , labelled diagrams [3l| and sequences of 
integers [27], [32|, |33| . However, these approaches could not be generalized to quasi-degenerate systems because they 
dealt with denominators (eo — ej) n with n > 1, whereas more complex denominators occur in the quasi-degenerate 
case. Moreover, they could not be used to discover efficient resummations because they do not capture the recursive 
structure of the RS series. 

To find out which combinatorial object is best suited to deal with the RS series, we had to use a rather elaborate 
algebraic analysis of the time-dependent perturbation theory [l4 j, 34] . It turned out in the end that planar binary trees 
provide the most faithful representation. As we shall see, the relation beween a tree and the corresponding term of 
RS series is straightforward and the recursive structure of the terms is transparent. The faithfulness is demonstrated 
by the easiness of the proofs and by the relation between geometrical properties of the trees and analytical properties 
of the terms, for instance for resummation. 



B. Binary trees 



We first give examples of the planar binary trees that we use. We denote by Y n the set of planar binary trees with 
n inner vertices. A tree is called binary if each vertex has either zero or two children. It is called planar if two trees 
are different when they can be deduced one from the other by moving one edge over another one. For example, the 

two planar trees | and } are different. 

There is a single tree with zero inner vertex: Y = { \ }. There is one tree with one inner vertex Y\ = { Y }, two 

trees with two inner vertices Y2 = { ^ , }, five trees with three inner vertices Y3 = { , ^y" , "Y* , Y^ , Y^ } 

and fourteen with four inner vertices Y i = {^^Y, ^Y,^T,^ ^^Y^j^^j^^iY^,^,^ ^ >1 >!}• The vertical line 
of a tree is called the root. The trees of Y n (with n > 0) can be built from smaller trees by the following relation 
Y n = {ti V t2 ■ ti S Yfc, ti G F„_fc_i, k = 0, . . . , n — 1}, where t\ V ti is the grafting of trees t\ and £2, by which the 
roots of t\ and ti are brought together and a new root is grown from their juncture. Pictorially: 
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For example, V|=Y,|vY=T7Yv'Y = T ■ Except in figure [TJ the vertices are not drawn explicitly 
for notational convenience. The inner vertices of a tree are the vertices to which three edges (or two edges and the 
root) are incident, its leaves are its vertices to which a single edge is incident. A leaf is oriented either to the left 
or to the right. Each tree of Y n has n inner vertices and n + 1 leaves. The order \t\ of a tree t is the number of its 
inner vertices. If C n denotes the number of elements of Y n , the recursive definition of Y n implies that Co = 1 and 
c» = E fc=o CfeC-»-fc-ij so that C n is a Catalan number. The recursive definition of planar binary trees make them 
very easy to generate with a computer. 

C. Recursive relation between trees and RS terms 
In ref . IbH . we showed that the wave operator can be written as the sum 

oo 

n=l teY n 

where 

(It = ]>>}^<j|, (4) 

with \j) an eigenstate of H in the model space and \i) out of it (i.e. P\j) = \j) and P\i) — 0). The scalars w t iJ are 
defined by an exceedingly simple recursive relation: for any tree t different from | , there are two trees t\ and t 2 such 
that t = t 1 Vt 2 . Then, 

%(k\V\l)' Ji 



y -^^t 2 (5) 



k,i e i ei 

with the special cases lu^ — —S^ if t\ = | and uj\^ = +8i.j if t-z = | . It is clear from the definition that (It is built 

from a product of \t\ matrix elements of V. Note that, when t\ = \ , both i and k in ajf* correspond to eigenstates 

outside the model space, while when £ 2 = I , both i and k in ujf^ correspond to eigenstates in the model space. 
We now prove that the expansion \ = St^| is a solution of eq. @. First notice that taking the commutator 

[(It, Ho] amounts to multiply uj 1 / by ej — ej. Let t = t\ V t%. Then, [(It, Ho] is equal to QVP, or QV(lt 2 or —d tl VP 
or — (l tl V(lt 2 , respectively, according to whether t\ = t 2 = |, or t\ = and t 2 ^ , or t\ ^ | and t 2 = |, or t± ^ | and 
t 2 =/= |, respectively. The result follows by summing over all trees t\ and t 2 . 

This very simple proof illustrates the fact that non trivial results can be easily obtained because trees encapsulate 
the recursive structure of the RS series. 

A more explicit version of the recursive relation ((5|) will sometimes be needed. For t = tiVt 2 , we define equivalently 
(l t by 

n, = E m ^fW\ if tl = | and fa = I , (6) 

< I f> J — p ■ 



'J 



-J 



n, - -Y. ^ ma " Pnm '< if^l a»dfa=[, (8) 

V J 

where the exponent of \ j p ) and l*^} means that state \ j) is a basis state of M and \i) a basis state of the complementary 
of M. 
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D. Examples 

As examples, we give all fit for \t\=2 and 3, the value of fit for |t| = 1 (i.e. t = Y ) being given by eq. ([5])) 

Qt = _ i^)(^i^Nr)fa p i^K3 p )(»n for t = ^ 5 

Qt = y tf)ti\V\$)f\VMtf\ for t=fj 
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fit 



y l^)^|y|^)(^l^l^)(^|y|0<^l for t= f_ 
. . . (ej 4 — ej 3 ) (ej 4 — e^ 2 ) (ej 4 — ) 



By looking at these examples, one may wonder whether some denominators could take the value zero. We show now 
by induction that this never happens. In other words, fit is never singular. This is true for t = Y because, in eq. (f6j). 
|i) belongs to M and \ j) does not belong to M. Thus, |e,; — ej\ is greater than the gap between the model space and 
the rest of the spectrum. Assume now that no fit is singular for trees with \t\ < n. Take a tree t with \t\ = n + 1. 
Then, t = t\ V £2, with \t\\ < n and fal < Assume that neither ii nor t% is equal to the root | . Then, eq. ((9]) and 
the previous remark about |e$ — e 3 -| shows that fit is not singular. The same is true if t\ or t% is equal to the root. 
Thus, the denominator fi t can never take the value 0. Again, the proof is made very simple by the recursive structure 
of f2 t . 



E. Direct relation between trees and RS terms 



The recursive relation ([5]) is very useful to derive proofs, but we also need a non-recursive expression for the terms 
of the RS series. We now give an explicit relation between t and fit [141 ], This construction can be followed in figure Q] 



for t = 1 . Consider a tree t with \t\ — n and number its leaves from 1 for the leftmost leaf to n + 1 for the rightmost 
one. The numerators in the expansion of fit are \ii)(ii\V\i2) ■ ■ ■ (i n \V\i n +i)(i n +x\. The state \ik) belongs to M (i.e. 
is i P }) if leaf k is oriented to the right and does not belong to M (i.e. is if leaf k is oriented to the left. In the 
example of fig. [TJ the denominator is | if } (i f | V | i% ) {i% \ V \ i f } (i f | V | if ) (i% \ V | if ) (if | 

Then, for each inner vertex v of t, take the subtree t v for which v is just above the root. In other words, t v is 
obtained by chopping the edge below v and considering the edge dangling from v as the root of t v . In the example of 
fig. [TJ we have four inner vertices labelled a, b, c, d. For vertex a, the subtree t a is the full tree t. For the other inner 
vertices the subtrees are proper (i.e. different from t) and are given in fig. [TJ For each tree t v , denote by l(v) the 
index of its leftmost leaf and by r(v) that of its rightmost leaf. Then, the denominator is the product of e< r .. — e^,, . , 
where v runs over the n inner vertices of the tree. In the example of fig. [TJ t a gives the denominator (ej 5 — e^), 
gives (ej 3 — e^), f c gives (ej 2 — 6^) and ^ gives (ej 5 — ej 4 ). Finally, the whole fraction is summed over . i n +i and 
multiplied by ( — where d is the number of leaves pointing to the right. It can be easily proved that the above 
construction satisfies the recursive equation ([5]) . If we take the example of the tree t of fig. [TJ and consider for example 
the inner vertex labelled by 6, we get a contribution (e^ 3 — e^) to the denominator. The total term is 



fit - (-I) 2 ^ 



\i?)(^\V\i?){ii\V\if)(if\V\i^\V\if)(if\ 



. . . . . {&i 5 e ii)( e ia e ii)( e i3 e ii)( e i5 e u) 

L l L 2 L Z L A L b 



FIG. 1: A tree with its labelled leaves and its three proper subtrees. The inner vertices are labelled a, b, c and d. The couples 
(l(v),r(v)) are (1,5), (1,3), (1,2) and (4,5) for v = a,b,c and d, respectively. 



As far as we know, this tree representation provides the first description of the general term of the RS series for 
quasidegenerate systems. 

III. RESUMMATIONS 

The tree representation is useful, not only to prove properties of the RS series, but also to derive new resummations 
of it. Indeed, any way to write the set of trees as a composition of subtrees gives rise to a resummation of the RS 
series. 

A. Summation over left combs 

For any tree t, we define the sequence t n by to — t, t n+ i = t n V |. For example, if t = , we have t\ = f , ta 
For any integer n, t n is obtained by grafting t on the leftmost leaf of a tree that has a single leaf oriented to the left 
and n leaves oriented to the right (such a tree is called a left comb whereas t n , n > is called a left comb grafting. 
Notice the uniqueness rewriting property: an arbitray tree u can be rewritten uniquely as t n , where t = | or t ^ v V |. 
The sum over graftings on left combs is made by defining 



n=l 

We can now calculate Vt' t in terms of Vl t . We first define 



3 J 



Then, we use Kvasnicka's trick and denote by Qi = \i^}{i^\ the projector onto the eigenspace of Hq with eigenvalue 
outside the model space, so that Q = Qi and eq. © can be rewritten f2 tl = — J^i Qi^tPVG^(ei). By repeating 
this argument we obtain 



n' t = ^QiSltfP '- PVG° P ( ei ) + PVG P {ei)PVG Q P {ei 

i 

= ^2 Qi&t(P + PVG P (ei)j ~\ 



The map P + PVG P (ei) goes from the model space to itself and the inverse is taken within the model space. The 
dimension of the model space is generally small and the numerical calculation of the inverse is fast. We transform 

(p + PVGpie^y 1 - ({H P - e t P + PVP)G P {e l )y 1 = (H Q P - e % P)(PHP - e.Py 1 , 
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where PHP = P(Hq + V)P and the inverse is again from M to M. To conclude, we assume that t = it V v, where u 
and v are different from | . Then, 

Q' t = ^QiSlvVSlvGPpieiXPHo - e.P^PHP - e^P)- 1 

i 

= y2 Qi n u vn v (PHP - e iP r = E mn _ uVnvim , 

where \j) and ej are the eigenstates and eigenvalues of PHP, so that P = |j)OI- In other words, summing over 
all left combs amounts to replacing the eigenstates of PHqP by the eigenstates of PHP. 



B. Accelerarated summation over left combs 



The formula we obtained for Q' t suggests to look for another expansion of the RS series involving only the eigenstates 
and eigenvalues of P(Hq + V)P. The uniqueness rewriting property implies that il = P + Q' + J^t ^'t> where the 

sum is over all trees t with two or more inner vertices and such that t ^ u V | , but this resummation can still be 
improved. 

We rewrite the Kvasnicka-Lindgren equation. 

\x,h ] = J2^ Ho - Ho J2^ = J2^ H °- e ^ 

i i i 

Thus, eq. © becomes £\ QiX(H ~ e<) = QVP + QV X ~ xVP - X V X that is, 

E QiXP(Ho ~ ei + V)P = QVP + QV X - X V X , (10) 

i 

or, for all the eigenvalues of Hq outside the model space: 

Q.X-P = QiVG P (ei) + Q l V X G P {e i ) - Q lX V X G P {e l ), 

where 

17) (71 



G P {z) = (PH P + PVP - zP)- 1 =E 



Cj — z 
3 J 

This equation is similar to the Kvasnicka-Lindgren equation and can be solved graphically and recursively by the 
same process. Notice that the leading term of the recursive expansion of X is now il' = J] QiVGp(ei). 

Y _ i 

Let us call a tree t right-normalized if there is no edge in t such that the tree t' obtained by pruning the tree t 
at that edge can be written t\ V | with t% ^ \ . The only left comb possibly contained in a right-normalized tree is 
Y ■ The number of right-normalized trees is considerably smaller than the number of trees. For n=l, 2, 3 and 4 the 

number of right-normalized trees is 1 ( Y )> 1 ( ), 2 ( , Y^ ) an d 4 ^Y^^S^'tY^" Y^) instead of 1, 2, 5, 14. The 
right- normalized trees are enumerated by the Motzkin numbers [35j | . 

The solution of eq. (|10l ) is then SI = P + f2f, where t runs over right-normalized trees and f2t can be defined 
recursively by 

tit = E ®f V ®® ifti=l and* 2 =|, (11) 

ft, . Y. *t>-\«it,+ \, (12) 

fk = ~ E mi ^ im iftx/l aadfe/l. (13) 

i£M,jeM 3 1 

Once again, this approach reduces the number. It can be expected that this resummation considerably accelerates 
the convergence of the series. Indeed, when the model space is well chosen, the matrix elements of QVP are smaller 
than those of PVP. In this expression, all terms involving powers of PVP have been resummed. 
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C. Alternative accelerated tree expansion 



A similar resummation, that we omit, can be obtained by summing over right combs. However, it is not so practical 
as the previous ones because it is usually not easy to invert QHQ — Qz. If this inversion is possible, then we can 
simultaneously sum over right and left combs, at least if we assume that a finite gap exists between the eigenvalues 
of QHQ and PHP. We make this assumption, proceed as before and transform the Kvasnicka-Lindgren equation. 
We define Pj — the projector onto the eigenspace of PHP with energy e,j and use similar notations for the 

projectors and energy levels of QHQ. Thus, P = Y^j Pj and we can rewrite the lhs of eq. (p~0|) as 

QiXP(Ho -e, + V)P = J2 QiXP(Ho - t* + V)Pj = £ QiX& ~ ei)P s = £ Q(e 3 - H Q ) X P 3 . 

i ij ij j 

Thus eq. (JTOj) becomes J2 3 Qfa - H )x~Pj = QVP + QV X ~ xVx, or 

J2Q(e j -H -V) X P j = QVP-xVx- 

j 

For all eigenvalues ej of PHP, this gives us 

xFi = Sfa)VT j -S& j ) x Vx'Pj, 

where we set: 

S(z) = (zQ-QHQ)-\ 

We can write directly the wave operator as a sum parametrized by trees: Q = P + ^2t>\ wnere is defined 
recursively for all trees by 

fit= E mvm\ = ^ s(JSj)QV p. tft=l> (14) 

i£M,j£M 3 1 jeM 

Sh = - X) m ^™f m = - E stefavntjPs ift = * 1 v* a . (is) 

As for the one involving a direct resummation over right combs, this resummation is not as practical as the previous 
one because the eigenvalues of QHQ are generally not known. However, this algorithm provides a way to express the 

(A C\ 

eigenvectors of a matrix M = ( ^ ^ j in terms of the eigenvalues and eigenvectors of A and B. 



D. Fourth resummation 



We use now another representation of trees in terms of trees. Since any tree t decomposes uniquely as t — t\ V ti, 
a recursion on the left hand side of the decomposition shows that there is a unique integer k and k unique trees Vi 
(with i = 1, . . . , k) such that 

t = Ck(v%,. . . ,v k ), 

where 

C k (v 1 ,...,v k ) = ((...(I Vt; 1 )V...)v« M )v»fe. 

In words, C k {v\, . . . , v k ) is obtained by taking a left comb with k right leaves and grafting v k on the lowest right 
leaf, v k -i on the next one, up to v\ on the highest right leaf. 
If we sum over k and over all trees Vi we find 



oo 

F = | +Y,£k{F,...,F), 
k=l 



(16) 
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where F stands for the formal sum of all trees. 

We can plug the expansion of the wave operator of section fill Cl in this equation and get for t = Cf.(vi, . . . , Vf.): 

n t = (-l) k £ S(e Jk+ J...S(e^QVP n Vn vi P n V...Vn Vk P ]k+1 . 

jl---jk + l 

Therefore, eq. (|16j) gives us the non-perturbative equation 

oo 

X = £Sfo)VP;+£(-l) fc £ S(e jk+1 )...S(e jl )QVP jl VxP h V...V X P jk+1 . (17) 
j fc=i ii— jfc+i 

If we use instead the first expansion Eq ([3]) of the wave operator and plug the corresponding recursive formulas in 
Ck{vi, . . . ,v k ), we get 

Qt = £ So(e Jk+1 )...S (en)QVPnVn vi P j2 Vn v2 ...Vn Vk P Jk+1 

jl—jk+l 

where fi := P and So(z) := (zQ — QifnQ) -1 , so that: 

oo 

n = P + £(-1)^ ^2(...((QiVGp(ei))VQ,Gj,(ei))...VQGp(ei)). 
fe=i i 

We believe that these relations are new. Similar expansions would follow by considering the symmetric expansion 
of trees (using the grafting operation of a tree on the right-most leaf of another tree instead of the grafting on the 
left-most tree). 



E. Relation with previous works 

The previous works on quasidegenerate perturbation theory correspond to a summation which is symmetric to our 
first alternative expansion (in the sense that they focus on the operators (a — QHQ)). By following a line suggested 
by Kvasnicka and Lindgren flU HH , several authors transformed the Kvasnicka-Lindgren equation into 

£(e, - QHQ) x Pi = QVP-xVP-xVx- 

i 

A resummation (similar, but slightly more involved than the one leading to eq. (fT7|)) gives the equation [36l |37| 

oo 

X = 5>( e ,)^P,+£(-l) fe £ S(e jk+1 ) ...SiejJVP^V + V X )P h (V + V X ) .-.(V + V X )P jk+1 , 
j k=l h—Oh+x 

which is a generalization of the degenerate case [3^, |39|. Suzuki and Okamoto further solved this for X t> u t their 
rather complex result was not applied to concrete problems, as far as we know. Up to a left/right symmetry, their 
result is similar to our first alternative expansion, excepted for the fact that PVP and PV S{e)V P are grouped in a 
single term. 

The main practical difference between the result obtained by Suzuki and Okamoto and our second alternative 
expansion is the fact that we do not consider the term PVP as a perturbation, we treat it exactly. This is important 
because, when the model space is well chosen, PVP is larger than the non-diagonal terms PVQ. 



IV. GREEN FUNCTION OF DEGENERATE SYSTEMS 



In this section, we discuss a question related to the Green function of degenerate systems. Consider a Hamiltonian 
H with a degenerate energy e . The eigenstates of H with energy en span a vector space M. The projector onto the 
model space M is denoted by P. In a series of recent papers 10 ! 1 J , we proved by non-perturbative methods that 
there are eigenstates \i) of Hq, called the parent states, such that the usual Gell-Mann and Low wavefunction has a 
well-defined limit when the adiabatic parameter e goes to zero: 

|w G ml) = hm 



>o (i|E7 £ (0,-oo)|i)' 
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where U e (t,t') is the evolution operator in the interaction picture. Morover, we showed that the parent states \i) are 
eigenstates of H (with energy e ) and are also eigenstates of PVP. As a consequence, (i\V\j) = for i ^ j if \i) and 
\j) are parent states. 

Since the parent states solve the problem in a non-perturbative approach, one might be tempted to use them in 
the perturbative one. In other words, we pick up a parent state, say |0), and we calculate the Rayleigh-Schrodinger 
series corresponding to the projector |0)(0|. However, as noticed by Toth [43|, a problem appears in the perturbative 

expansion. This problem can be illustrated by a simple example: for t = , we have [l4j 

Qf = y 

rrt (eo - &i + ie)(eo - e% + 2ie) ' 

where e is the adiabatic switching operator, which tends to zero at the end of the calculation. Now, if \i) belongs to 
the model space, then (i|V|0) = (because i ^ 0) and the expression converges although e; = eo could have brought 
a problem. In other words, using a basis of parent states for M has made the expression convergent. However, this 

trick does not always work. Indeed, for t = f,we have [l4[ 

n t = y mvm\v\o)(Q\ 

* ijojito ( e ° " e J + * £ ) ( e o - ei + 2ie) ' 

If \i) belongs to the model space and \j) is out of it, then we have e^ = eo and the limit e — > is not defined 
because nothing insures that (i\V\j) = or (j\V\0) = 0. In other words, the convergence problem was solved at the 
non-perturbative level but remains at the perturbative one, so that the perturbative expansion has to be resummed 
in a proper way. 

Now we show how to solve this problem by using a trick related to the Hamiltonian shift proposed by Silverstone [3] ■ 
We assume that e\ — + are nondegenerate. Then, we rewrite H = Hq + V = Hq + V' , where 

H' Q = H + J2\i)(i\V\i)(i\, 

We build the RS series for Hq and V with the one-dimensional model space M spanned by |0). Thus, P' = |0)(0|. 
This gives us P'V'P' = 0. As a consequence, Ct t = if t = t\ V | with t\ ^ | . We can write Q' = Qo + Q, where 
Q is the projector corresponding to the initial problem and Qo = P — P' is the projector onto the basis states of 
M different from |0). Then, we have QqV'P' = and QV = QV. This gives us Q'V'P 1 = QV'P' = QVP', which 

simplifies the evaluation of f} { for t = Y • Similarly, P'V'Q' = P'VQ simplifies the evaluation of Q t for t = t\ V *2- 
Finally, the identity Q V'Q = gives us Q'V'Q' = Q VQ + QVQ + QVQ for the evaluation of Sl f with t = \ V t 2 . 
This gives us the following recursive expression for t — t\ V ti- 



= E i^^iy(Q + Q )a 2 io)(o| +E i^)(^ga 2 |o)(o| . {ti=] 

fl t = if *i ^ I and t 2 = 1 , 

n( = y \i Q )(i Q \n tl \o){o\VQn t2 \o){o\ _ y \iQ°){iQo\n tl \o){o\vQn ta \o){o\ . f ^ ^ t 

i e ej . c e i 

The recursive expression shows that all terms are well defined if all e[ — + (i|V|i) are different and if e is different 
from the energies ej out of the model space. 



V. CONTINUED FRACTIONS 



We discuss here some continued-fraction resummation of the RS series. Such a (generalized) continued fraction 
formula was found to be very efficient for calculating nuclear properties [45| . The combinatorial structure of continued 
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fractions was studied in detail by Flajolet [46|, |47[. On the other hand, Lee and Suzuki derived a continued fraction 
expression for \ and the effective Hamiltonian for a degenerate system [48l . l49j . Other implementations of continued 
fractions for pertubation theory can be found in the literature (50 - 4541 ) . 



A. The Suzuki-Lee formula 



Suzuki and Lee [49l ] start from the Kvasnicka-Lindgren equation 

[ X ,H ] = QVP + QVx-xVP-xVx- 

For a degenerate system PH = eoP, so that xHo — eoX an d (eo — H )x = QVP+QVx~ X^P ~xVx- This equation 
is then reordered into (eo — QHQ + xVQ)x = QVP — xVP- They consider the iterative equation (see eq. (3.27) 
p. 2102). 

{e -QHQ + Xn-iVQ)xn = QVP- X n-iVP, 
with the boundary condition xo — 0. In other words 

Xn = (e -QHQ + X n-iPVQ)-\QVP-Xn-iVP). 



B. A new continued fraction expansion 

The Suzuki-Lee formula has two drawbacks: it is restricted to degenerate systems (55| and it requires the inversion 
of eo — QHQ + Xn-iPVQ, which is usually infinite dimensional. To solve these two problems, we start from eq. (fTO]) 
and, by using QiQ — Qi, we rewrite it QtxP(H — a)P — QiVP + QNx ~ QiX^X- We transform this equation into 

Q iX (P(H- ei )P + PV X ) = QiVP + QiVx- (18) 
Thus, we define the system of recursive equations 

Q l Xn = {Q l VP + Q l Vx n -i){PHP-e l P + PVxn-ir\ 

Xn ^ QiXm 

i 

with the boundary condition xo = 0. This generalized continued fraction has convergence properties similar to that 
of Lee and Suzuki, it is well-defined for quasi-degenerate systems and the inverse is computationally easier because it 
is done within the model space. 



VI. BIJECTIONS 



Many other combinatorial objects have been used to represent the terms of the RS series in non-degenerate or 
degenerate cases. Each of these representations is useful for specific applications. It is therefore important to describe 
the relation between the most important of them (Bloch sequences, Dyck paths, braketings and non-crossing partitions) 
and the trees. Most of these representations are valid for degenerate systems. Thus, we start by presenting the first 
terms of the RS series of degenerate systems, where eo is the energy of the states of the model space. We also give an 
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TABLE I: Correspondence between several representations of the RS terms for n=l,2,3: numbered tree, Bloch sequence, Dyck 
path, bracketing, non-crossing partition and item in Olszewski's list of examples [281 ]. 



operator representation of these terms, with R = Q(eo — Hq) 1 Q. 

at = ^\mmm =RVP for «_ Y , 

e - e it 
•f-f (e - e i2 )(e - e^j 



fi, 



n* = 



y \§Mm^M\MM^Ml Rzypypyp for *=> 

( e o-eu)(eo-e I1 )(e -e n ) 

_ E i^^i^^imy = _ R2 y RVPVP for t= ^ 

e - e 8l ) (e - e l2 (e - e;J 

v i*?> </? i y-i»r > <«riv-i«?> i vi»r> <»r i = _ R * VPVRVP for t= y 



(e - e ll )(e - e ll )(e - e is ) 

"± «a 

n t = - V l'?)<'?l^l'?)<'?l^l^<»3nvi»f Kifl = _ RVR > VPVP for t= f 

. ff. e - e n ) e - e l2 e - e i2 ) 

^ = E \§m^m^mm = RVRVRVP for t= 

. j-f. (e - e l3 ) (e - e i2 ) (e - e u ) 



A. Bloch sequences 

Bloch [27j was the first to write the general term of for degenerate systems. To describe his result, we consider 
the wave operator fi(A) of the Hamiltonian H + XV. We have the series expansion 

00 

n{\) = p + ^2x n n n , 
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TABLE II: Correspondence between several representations of the RS terms for n=4: numbered tree, Bloch sequence, Dyck 
path, bracketing, non-crossing partition and item in Olszewski's list of examples [28l ]. 



with 

n n = S (kl) VS {k2) V ...VS (kn) VP, (19) 

k\ .. . .,k n 

where S° — —P, — R k for k > 0, and where ki,...,k n run over the Bloch sequences. A Bloch sequence is an 

n-tuple (fcife . . . fc„) of non-negative integers, such that fei + • • • + k m > m for m < n and k\ H h k n = n. The Bloch 

sequences for n=l, 2 and 3 are given in table U and for n=4 in table [TIJ These combinatorial objects are enumerated 
by Catalan numbers (see Example 6.24 p. 1 80 of ref. [25[ and item m 5 in Stanley's Catalan addendum). Since Bloch's 
publication 27 1, they were widely used |5a - l59| to represent the general term of the RS series for degenerate systems. 



The bijection between trees and sequences is defined by 4>( Y ) := (1), </>(s V | ) :— K(j>(s), (/>(] V s) :— (1) • cj)(s) 
and (f>(s V t) := Kcf>(s) ■ <f>(t) where s ^ | , t ^ | , where the product on the right hand sides is the concatenation 
of sequences (i.e. (k\ . . . fc m )(fc m +i . . . k n ) = (fci . . . fc n )), and where the operation K acts by K{k\, k 2 , ■ • ■ , k n ) = 
(ki + 1, &2, . . . , k n , 0). Notice that </> is clearly one to one (injective). The fact that it is actually a bijection follows 
e.g. from the fact that trees and Bloch sequences are both enumerated by Catalan numbers. 

We prove now its compatibility with the tree and Bloch expansions by induction. We write t = t\ V ti and we 
consider the usual four cases. The first case is t — Y , ti = t 2 = \ ■ This is the starting point of the inductive proof. 
We have fit = RVP = S 1 VP = f2i, where we used the fact that the only Bloch sequence for n = 1 is (1). Let 
us denote by Qb the contribution of Cl n in eq. (|19|) corresponding to the Bloch sequence B = (ki ...k n ). Thus, we 
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showed that, for t = Y an d B = (1), we have Vl t — Qb- 

Now, assume that fl t = for all trees t with \t\ < n and choose a tree t such that \t\ = n + 1. Then, t = t± V ti 

with \t\\ < n and 1^1 < We consider the three remaining cases. If t\ ^ | and t% — | , then f2 t is obtained from 
eqs. Q and ([5} 51 t = —RD, tl PV P . By the induction hypothesis, we have 57 tl = Sl^ for B\ = 4>(ti). Therefore, 
fl t = RQ, Bl (-P)VP. If Bx = (ki ... k n ), then O t = ft B with 5 = {k\ + 1, fc 2 , . ■ • , fc n ,0) = KB\ and the property is 

proved. The second case is t\ ^ | and f 2 ^ | . Then eqs. ((U) and (HJ) give us Q t = R^t 1 (^P)VQi} t2 - If </>(*i) = 
B x = (ki . . . k p ) and 4>{t 2 ) = B 2 = {h... l q ), then Q t = B with 5 = (fcj + 1, fc 2 , . . . , k p , 0,h,..., l q ) = K(B 1 )B 2 . 

Thus, (j>{ti V t 2 ) = K(B 1 )B 2 = K<t>{l\) ■ (j)(l 2 ). The last case is h = | and t 2 ^ \ . Equations g]) and © give us 
fit = RVQtt t2 . If B 2 = <f)(t 2 ) = (h... k n ), then Q t = Q B with B = {lk x . . . k n ) = (1) • <j)(t 2 ). 

B. Dyck paths 

In ref. [27|, Bloch also defined geometrical objects that are Dyck paths rotated by n/4. We prefer to use Dyck 
paths because they are thoroughly studied in the combinatorial literature (see, for example Item i, p. 221 of ref. |26j 
or Example 6.2, p. 151 of ref [25j, where they are called mountain ranges). 

The bijection between Dyck paths and Bloch sequences used in ref. [271 ] is well known in the combinatorial litera- 
ture [111 p. 168 and 181]. To build the Dyck path corresponding to the the Bloch sequence (kik 2 . . . fc„), start from 
the origin and make k± steps in the North-East direction, then make one step in the South-East direction, then k 2 
steps in the North-East direction, then one step in the South-East direction, and so on. This bijection is illustrated 
m tables U and |TT] 

C. Non-crossing partitions 

The terms of the RS series of degenerate systems can also be described by non-crossing partitions. This corre- 
spondence was studied by Olszewski [28[ because it leads to useful factorizations of the RS terms for non-degenerate 
systems. 

Consider the examples at the beginning of section IVI1 If we just look at the denominators, we see that each of 
them can be deduced from the right comb by saying that some indices are equal. For the five trees with \t\ = 3, the 

term corresponding to t = | / can be obtained from that of t = \ by stating that = e !2 = . The other terms 
follow from = e,- 3 , = ei 2 and ei 2 = ej 3 . 

More generally, for a given tree t, we say that two indices j and k are equivalent if and only if ej. = in 
the denominator of £l f . The sets of equivalent indices form a partition of {1, . . . ,n}. We recall that a partition 
Bi\B 2 \ . . . \Bk of {1, . . . , n} is a set of disjoint subsets Bi of {1, ... , n} whose union is {1, . . . , n}. Each subset Bi is 
called a block of the partition. The partitions corresponding to the RS terms of order 1 to 3 are given in table U and 
for order 4 in table [TTJ It will be shown that these partitions are non-crossing. Two blocks A and B of a partition are 
said to be crossing if there are a, b in A and x, y in B such that a < x <b < y or x < a < y < b. A partition is called 
non-crossing if no two of its blocks cross. 

We build by induction a partition from a tree. Note first that for tree t of order n we deal with partitions of 
{1, . . . , n}. A tree of order n has n + 1 leaves. Thus, the index of the rightmost leaf is not used in the partition. 

For a given tree t = t\ V t2, we call P (Pi, P 2 , respectively) the partition corresponding to t (ti, t 2 , respectively). 
We consider the usual four cases, (i) For t\ — t 2 = \ we associate the partition P = |1|. (ii) For t\ ^ \ and t 2 = \ , 
eq. ([8]) gives us the additional denominator eo — e±. Therefore, an additional index is equivalent to 1. This index is 
that of the rightmost leaf of ti, which was not used in P x . Therefore, the partition P of t is obtained from P x by 
adding to the block of P\ containing 1 the index of the rightmost leaf of ti, which is |ii| + 1. (iii) For t\ — \ and 

t% 7^ | , eq. (O gives us the additional denominator eo — e\. However, the leaf denoted by 1 is new and no block of 
P 2 should contain it. Therefore, the partition P of t is obtained from P 2 by increasing all numbers of P 2 by 1 and by 
adding the block |1|. (iv) For t\ ^ | and t 2 ^ | , we compose the previous cases. We build P by adding the index 
\ti \ + 1 to the block of Pi containing the index 1 and we increase all indices of P 2 by \ti \ + 1. We check that, in all 
cases, the resulting partition is non-crossing. 

Note that Olszewski [28| does not explicitly use non-crossing partitions. He describes the general term of the RS 
series for non-degenerate systems by drawing a circle with n points and pinching some of these points. The relation 
with non-crossing partitions is straightforward: all points that are pinched together belong to the same block of the 
partition. However, his correspondence between partitions and terms of the RS series is not the same as ours. 
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D. Bracketing 

Following a suggestion by Brueckner (60[, Huby and Tong [29],[30[ proposed, for the energy of nondegenerate systems, 
a solution in terms of bracketing, that is isomorphic with Stanley's problem e 5 in his "Catalan addendum" . 

Consider the example of t = . We have Qt — R 2 VPVP. For a nondegenerate system P = |0)(0| and f^rlO) = 
R 2 V\Q){Q\V\Q). The rule of the game is now to insert expectation values to disjoin powers of R. Thus, Qt|0) = 
—R(0\V\0)RV\0). We use Tong's pictorial representation, where R is replaced by *, V by o, (0| by ( and |0) by ), so 
that f2x|0) = — * (o) * 6). The bracketings for |t| = l,2 and 3 are given in tableland for order 4 in table[llj Brueckner's 
bracketing is a powerful way to simplify the RS series, for instance by including the vacuum expectation value (0| V"|0) 
into Ho, so that all RS terms involving it cancel. More details and examples can be found in some textbooks [61| or 
review papers [62~| . 

We now describe the connection between trees and bracketings. Denote by b(t) the bracketing corresponding to 
t in Tong's representation. Assume that b(t) is known for all trees t with \t\ < n and take a tree t = t± V of 

degree n. It ti — | and ta ^ I > then f2t|0) — i?VT2t 2 |0), so that b(t) = *ob(t2). It t\ ^ | and t 2 = J , then 
n t |0) = -Rft tl \0){0\V\0) = -R(0\V\0)n tl \0), so that b{t) = -*{o)b(h). It ti ^ | and t 2 j= \ , then Q 4 |0) = 
-Rn tl \0)(0\Vtt t2 \0) = -R(0\Vn t2 \0)n tl \0), so that b(t) = -*{ob(t 2 )b(ti). Note that this bijection is different from 
the one used by Tong [30| . 

An equivalent bijection is obtained by numbering the inner vertices of the trees. The operation v that associates 
to each tree t its numbered tree is defined as follows. We denote by i/(t)[fc] the numbered tree obtained from v(t) by 
adding k to all the vertex numbers. Then, v(t) can be defined recursively. If t = , then v(t) — \ (no number). If 
t = Y i then v(t) is obtained by assigning the number 1 to the root. If t = t\ V t2, then the root has number 1, the 
inner vertices of ti (if any) are numbered as i/(f2)[l] and the inner vertices of t\ (if any) are numbered as ^(ti ) [|*2 1 + !]• 
The numbered trees for \t\ < 3 are given in table table U and for \t\ = 4 in table ITIl For any tree t, we can build the 

V . 

sets of numbers belonging to the same line oriented to the left. For example, if t = J , the numbering isj 1 and the 
sets are |1|24|3|. These sets form a non-crossing partition which is the same as that used by Olszewski [2g|. To obtain 
b(t), first number the \t\ stars of b(t) from 1 to \t\ from the left to the right. Then, consider a block B = ki . . . k p 
of the partition. If p = 1, do nothing, if p > 1, then write a ( after star number k\, a ) before star number k p and 
replace star number ki with 1 < i < p (if any) by ) * (. 



VII. CONCLUSION 



Combinatorial physics is an emerging field that uses modern tools of algebraic combinatorics to solve physical 
problems. It was born with the investigation of the algebraic structure of renormalization in quantum field theory [63[ 
and showed its ability to deal with many-body problems [64|. 

We showed that combinatorial physics is able to solve such long-standing problems as time-independent perturbation 
theory. The RS series is at the heart of many applications of quantum mechanics. It is also equivalent to more 
sophisticated methods such as Feynman diagrams 65] . It is even related to Wilson's renormalization group (66j j . 

Our combinatorial methods provided easy resummations of the RS series. It remains now to test their convergence 
properties. 

Note that Arnol'd also used trees in perturbation theory [67[ . However, his trees are essentially different from ours 
because they describe the successive degeneracy splitting due to higher order terms [6o |. 
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